Resonant photon creation in a three dimensional oscillating cavity 
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We analyze the problem of photon creation inside a perfectly conducting, rectangular, three 
dimensional cavity with one oscillating wall. For some particular values of the frequency of the 
oscillations the system is resonant. We solve the field equation using multiple scale analysis and 
show that the total number of photons inside the cavity grows exponentially in time. This is also 
the case for slightly off-resonance situations. Although the spectrum of a cavity is in general non 
equidistant, we show that the modes of the electromagnetic field can be coupled, and that the rate 
of photon creation strongly depends on this coupling. We also analyze the thermal enhancement of 
' the photon creation. 
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I. INTRODUCTION 



< 

The existence of an attractive force between two perfectly conducting plates has been predicted by Casimir in 1948 
■ JjJ . It has been measured with accurate precision in the last years by Lamoreaux Q and Mohideen et al. || . These 
^ \ experiments confirm the existence of vacuum field fluctuations in the framework of field quantization with static 
. boundaries, and increase the interest in the case of dynamical boundaries as well. 

The dynamical effect consists in the generation of photons due to the instability of the vacuum state of the electro- 
magnetic field in the presence of time-dependent boundaries. In the literature it is referred to as dynamical Casimir 
effect [Q or motion-induced radiation || . Up to now no concrete experiment has been carried to confirm this photon 
generation, but an experimental verification is not out of reach. From the theoretical point of view it is widely ac- 
cepted that the most favorable configuration in order to observe the phenomenon is a vibrating cavity in which it is 
""^ \ possible to produce resonant effects between the mechanical and field oscillations. 

Many previous papers have focused their attention in the field quantization within a one dimensional cavity with one 
or two walls performing small amplitude oscillations, at twice the eigenfrequency of some unperturbed electromagnetic 
mode. For these cavities there exists a strong intermode interaction, which is a consequence of the equidistant character 
of the frequencies spectrum. The main features of the one dimensional model are that photons are created in all 
electromagnetic modes (due to mode- mode coupling), and that the total energy inside the cavity grows exponentially 
at the expense of the energy given to the system to keep the wall moving. A simple approach is to make a perturbative 
expansion in terms of the amplitude of oscillations, as was done in However this perturbative treatment breaks 
down after a short period of time due to the appearence of secular terms. In j?j the renormalization group technique 
was used in order to obtain a solution valid for a period of time longer than that of the pertubative case. There it 
was shown that the energy spectrum develops a non trivial structure formed by peaks travelling at the speed of light 
and bouncing against the walls (in agreement with other authors ||). Reference B makes use of the fact that two 
different time scales characterize the problem; the usual one, related to the wall's oscillation period, and a 'slow' one 
which accounts for the cumulative resonance effect. It is then posible to isolate the resonant part after averaging over 
fast oscillations the initial equations for the electromagnetic field modes. 

There are some works in the literature dealing with higher dimensional cavities . In fl(i[| the radiation emitted in 
each polarization of the electromagnetic field was computed perturbatively, when two parallel, plane, and perfectly 
conducting plates oscillate along the direction perpendicular to their surfaces. Such geometry constitutes the simplest 
example of an open three dimensional cavity. In pd| the authors obtained the distribution of the created photons for 
the case of parametric resonance inside a three dimensional cavity. In both cases jHJjyJ a perturbative method was 
applied so the results are valid in the short time limit. In B a nonperturbative analysis was presented, generalizing 
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the method of averaging over fast oscillations to higher-dimensional cavities. However, the intermode coupling was 
neglected, reducing the problem to that of one single parametric oscillator. 

Of particular interest is to find out how the finite temperature affects the photon production. This was studied in 
|l2f with a nonperturbative method (see also p3| ). A remarkable enhancement of the pure vacuum effect was found, 
but again neglecting the coupling between modes. 

In this paper we present a detailed analysis of the photon generation inside a three dimensional resonant cavity. We 
also discuss the finite temperature case, showing the enhancement of the effect with respect to T — 0. We apply the 
Multiple Scale Analysis (MSA) which provides us with a simple technique equivalent to summing the most secular 
terms to all orders in the perturbative treatment. In this way we can get a solution valid for a period of time longer 
than that of the perturbative case. We pay particular attention to the resonant coupling between different modes. 

The paper is organized as follows. In section II we obtain the time evolution of the quantized field by expanding it 
over the 'instantaneous basis '. For simplicity we deal with a scalar bosonic field. Following the steps given in fllij we 
arrive at an infinite set of coupled differential equations for the coefficients of the expansion. We also explain there 
how to compute the number of motion-induced photons for the zero temperature case. In section III we describe 
and apply the MSA to our problem. We find the coupling conditions between different modes that can be satisfied 
depending on the cavity's spectrum. In section IV we present a general analysis of the coupling conditions, and discuss 
some examples. In particular we find that the fundamental mode of a cubic cavity is coupled to another mode in 
the parametric resonance case, giving as result that the number of photons with two different frequencies increases 
exponentially in time. However, the production rate for the fundamental mode is only one half of that expected if we 
had neglected the coupling, as previous works did. At the end of this section we study slightly off-resonance situations. 
In section V we obtain an expression for the number of photons in each mode assuming that the field was initially in 
thermal equilibrium. Section VI contains our final remarks and comments on the generalization to the more realistic 
case of the electromagnetic field. 

II. SCALAR FIELD QUANTIZATION WITH MOVING BOUNDARIES 

We consider a rectangular cavity formed by perfectly reflecting walls with dimensions L Xl L y , and L z . The wall 
placed at x = L x is at rest for t < and begins to move following a given trayectory, L x (t), at t = 0. Note that we 
assume this trayectory as prescribed for the problem (not a dynamical variable) and that it works as a time-dependent 
boundary condition for the field. The field 0(x, t) satisfies the wave equation □</> = in 3 + 1 dimensions, and the 
boundary conditions 0| wa iis = for all times. The Fourier expansion of the field for an arbitrary moment of time, in 
terms of creation and annihilation operators, can be written as 

<Xx, f) = ^«;> n (x, f) +H.c, (1) 

n 

where the mode functions M n (x, t) form a complete orthonormal Q set of solutions of the wave equation with vanishing 
boundary conditions. 

When t < (static cavity) each field mode is determined by three positive integers n x ,n y and n z . Namely 



u n (x, t < 0) = J-f-Bin , -^sin ^y) J^sm ^z e^\ (2) 





(3) 

with the shorthand n = (n x , n y , n z ) 0. 

When t > the boundary condition on the moving wall becomes <p(x = L x (t), y, z, t) = 0. In order to satisfy it we 
expand the mode functions in Eq.(|l|) with respect to an instantaneous basis 



1 The inner product is the usual for Klein-Gordon equation, namely = — i /" cavit d 3 x [V'C* — ifiC] 



2 We are using units h ■ 
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„„<*,,> 0] = Eflfw^- (^) ^- ^ (If 2 ) (4 » 

= £<JL n) (<)*>k(*.M<)). (5) 



with the initial conditions 



In this way we ensure that, as long as L x (t) and L x (t) are continuous at t = 0, each field mode and its time derivate 
are also continuous functions. The expansion in Eq.(^|) for the field modes must be a solution of the wave equation. 
Taking into account that the <fkS form a complete and orthonormal set and that they depend on t only through 
L x {t), we obtain a set of coupled equations for Q^(t) fHf] : 

Ql n) +"l(t)Qi n] = 2A(;)5> kj Qj n) + A(t)]Tg kj Qj n) + A 2 (t) 5> k <?y Qj n) , (7) 

j j j,i 

where 



The coefficients <7kj are defined by 



and read 



L x (t) 



.9kj = L x {t) j dx -Qj-Vh ( 9 ) 



5kj = -aik = i ; ' Bjy z3z ., , io) 

10 if fee = J x . 

Furthermore, in deriving Eq.(^) we have used that gkiSji = L x J dx (dif] SL /dL x )(dLpj/dL x ), which follows from the 
completeness relation of the ip^s. 

The annihilation and creation operators a k and d k correspond to the particle notion in the 'in' region (t < 0). If 
the wall stops for t > tgnai, we can define a new set of operators, a k ut and d k ut , associated to the particle notion in 
the 'out' region (t > tfinai)- These two sets of operators are connected by means of the Bogoliubov transformation 

n 

The coefficients a n k and /3 n k can be obtained as follows. When the wall returns to its initial position the right hand 
side in Eq.(0) vanishes and the solution reads 

Q^\t > U = A^e^ + B^e-***, (12) 

with and B k being some constant coefficients to be determined by the continuity conditions at t = ifi na i. 

Inserting Eq.(|l2|) into Eqs.(|l|) and (@) we obtain an expansion of <ft in terms of d k and d k ' n for t > i nna i- Comparing 
this with the equivalent expansion in terms of a k ut and d k ut it is easy to see that 

a nk = V2^k£ k n) , PnU = V%^A£ ) . (13) 

The amount of photons created in the mode k is the average value of the number operator d k ° ut a k ut with respect 



to the initial vacuum state (defined through d k |0 in ) = 0). With the help of Eq.(ll) and Eq.(|l3[) we get 

(mo - (o,„ 1 4° ut dr i o in > = ^2^ k i4 n) i 2 . (w) 
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III. MULTIPLE SCALE ANALYSIS 



Up to this point the equations are valid for an arbitrary motion of the wall (we only assume L(0) = Lq and L(0) = 
because the wall is at rest for t < 0). We are interested in the number of photons created inside the cavity, so it is 
natural to look for harmonic oscillations of the wall which could enhance that number by means of resonance effects 
for some specific external frequencies. So we study the following trayectory 



L(t) =L (l + e$in(nt)+ef(t)), 



(15) 



where f(t) is some decaying function that allows us to meet the continuity conditions at t = (for example f(t) = 
—Slte~ at ). For small amplitudes of oscillations (e <C 1), the equations for the modes Eq.(Q) take the form 



L 



efl 2 sin(Ot) ^ ffkjQj n) + 2efl cos(ftt) ^ g^Q ( " } + +eO(f) + <3(e 2 ), 



(16) 

where O(f) denotes terms proportional to /, / and /. 

It is known that a naive perturbative solution of these equations in powers of the displacement e breaks down after 
a short amount of time, of order (efi) -1 . This happens for those particular values of the external frequency such 
that there is a resonant coupling with the eigenfrequencics of the static cavity. In this situation, to find a solution 
valid for longer times (of order e" 2 ^" 1 ) we use the multiple scale analysis (MSA) technique pl|. We introduce a 
second timescale r = et and expand as follows (we shall content ourselves with first order MSA) 



W(t) = Q^°\t, t) + eQ^\t,r) + 0(e 2 ). 



(17) 



The derivatives with respect to the timescale t read 



+ 3 2 Q< n)(1) ] 



(18) 



The initial conditions are 



j(n)(0) 



(0) = 



= <5n.k, 



Q(n)(0) 



(0) 



'\/— n ,k- 



(19) 



To zeroth order in e we get the equation of an harmonic oscillator 



W(0) ^[ n) (r)e' utt 



^k - ^k 
and using the initial conditions it follows that 



(20) 



4 n) (r = 0) = 0, 
b£\t = 0) = ^ 



-Sn 



k- 



(21) 
(22) 



To first order in e we obtain 



3 It is not straightforward to compute the next order corrections using MSA. The introduction of new time scales like n = 
et,T2 — e 2 t, etc, is in general not sufficient to determine the second order solution unambiguosly |L5|]. The renormalization 
group method |16fl (which is equivalent to MSA to first order in e) seems to be more appropiate to improve systematically the 
result. In any case, the next order corrections will be very small for e 2 t < Q _1 . 



4 



= -2%Q 



)(o) 



^) sin(m)Q k r 

,(n)(0) 



)(0) 



-ft 2 sin(Ot) ^ SkjQj )(0) + 2ft cos(ftt) ^ <?kj<9 t Qj 



(n)(0) 



o(/) 



(23) 



The basic idea of MSA is to impose that any term on the right-hand-side of the previous equation with a time 
dependency of the form e ±ZUJkt must vanish. If not, these terms would be in resonance with the left-hand-side term and 
secularities would appear. The terms contained in 0(f) are not relevant because they are exponentially suppressed, 
and do not produce secularities. After imposing that no term e +l " k * appear, we get 



dA 



(n) 



2 1,2 



(It 



2u> k L% 



ft. 



ft 



4 n) j(2 Wk - n) + + j) s (-^ -<»i + n)— g kj B 



(n) 



( w j + -u 3 -Q) + (wj - — )5(wk - cjj + ft) 



2wk' 

ft 
2w k 



(n) 



(24) 



In a similar fashion, the fact that no secularities should arise from the e %UJkt term leads to 



dB- 



(n) 



dr 



2 1.2 



TT 2 k 



f^A^6(2u; k ft) + £(-a/, + - + n)^-5kjA] 



ft 



2w k L 

E 



(n) 



( w j + - wj - ft) + (wj - ^)<5(wk - + 0) 



2w k * 

ft 
2w k 



(n) 



(25) 



The previous set of two equations are non trivial (i.e., lead to resonant behavior) if ft = 2uj k (resonant condition) 
Moreover, there is intermode coupling between modes j and k if any of the following conditions is satisfied 



ft 



There is an alternative, equivalent way of deriving the equations of motion 
that the solution of Eq. (n6|) is of the form 



QW(t) = 4 n) (*)e <Wkl 



flW(t) 



— iuJi^t 



(26) 
(27) 
(28) 

For e <C 1, it is natural to assume 



(29) 



where the functions and B k n ^ are slowly varying. In order to obtain differential equations for them, we insert 



this ansatz into Eq.(|l6|) and neglect second derivatives of A-j?' and £>! n) . After multiplying the equation by e ±l " k * 



,(n) 



we average over the fast oscillations. The resulting equations coincide with Eqs.(|24|) and (|25|). 

We derived the equations for 3+1 dimensions. It is very easy to obtain the corresponding ones in 1+1 and 2+1 
dimensions. In all cases the resonant conditions are given by Eqs.(^6|)-(^8|) above. The main difference between the 
1+1 case and higher dimensions is that in 1+1 the eigenfrequencies u>k are proportional to integers. The spectrum 
is equidistant and therefore an infinite set of modes may be coupled. For example, when the external frequency is 
ft = 2u>i, the mode k is coupled with the modes fc± 2. This has been extensively studied in the literature [j|j, , [L7| . 
In what follows we will be concerned with cavities with non equidistant spectrum. 



IV. RESONANT PHOTON CREATION 



In this section we shall solve the coupled Eqs.(24) and (|25|). We will see that there are different kinds of solutions 
depending both on the wall's frequency and the spe ctrum of the static cavity. Note that the spectrum is related to the 
cavity's dimensions through Eq.(^). In subsection [V A we will present a gen eral analysis of the resonant conditions 
and the solutions. We will show some particular examples in subsection IV B. 



5 



A. General analysis 



Let us consider the 'parametric resonance case', in which the frequency of the wall is twice the frequency of some 
unperturbed mode, say Q = 2cjk- Under this condition we expect that the number of created photons in the mode 
k will grow exponentially in time due to resonance effects. In order to find At* 1 ' 
we have to analyze whether the coupling conditions |u;k ± Wj| = 
resonant mode k will be coupled to some other mode j only if wj 
satisfied depending on the spectrum of the particular cavity under consideration. 



and from Eq.fl24|) and Eq.(|2jf) 

fl can be satisfied or not. If we set f2 = 2wk, the 
Wk = O = 2wk- Clearly, the latter relation will be 



First, let us assume that this condition is not fullfillcd. In this case, the equations for A± and B^ 1 ' reduce to 



dA 



(n) 



(It 



dB 



dr 



(n) 



— 1 / 7rfc 3 
2^k V Ac 

— 1 / 7rfc a x 
2wk V Ac 



B 



k ' 



.4 



(n) 



The solution that satisfies the initial conditions (|T|) and (|2^) reads 



(30) 
(31) 



B; 



(n) 



A 



k 

(n) 



: ^ k , n cosh(7A: :c r), 
sinh(7fc s r), 



(32) 
(33) 



where 7 = (k x /n)(ir/L x ) 2 . With the help of EqQ we obtain 

(TVk) = sinh 2 (7fc x T/), 



(34) 



where Tf — etf. In this uncoupled resonance case the average number of created photons in the mode k increases 
exponentially in time with a rate given by 2"/k x . The same result has been obtained in previous papers (cs. Ref. 
H and Ref. jjJI). There it was assumed that the coupling condition Wj = 3wk cannot be fullfiled for two and three 
dimensional cavities, essentially due to the non-equidistant character of the spectrum. As we shall see, this is not 
always true. In what follows we will solve Eq.(p4|) and Eq.(^5|) with coupled modes, and we will show some explicit 
examples. 

Let us now assume the existence of one mode, say j, in the infinite sum in Eq.(p3) and Eq.(p5|), which satisfies 



3wk- We obtain for A^ and B 



,(n) 



dA 



(n) 



dr 



dB 



(n) 



dr 



l(-k x B^ + (-iy*+ k * Jx A^), 



(-k x A { ^ + {-iy* +k * ]x B\ 



(n) 



(35) 
(36) 



where we have used that the relation = 3wk is equivalent to 



j 2 x = 9kl + 8 



L u 



(37) 



because the coupling coefficient <?kj is proportional to 8u y j v Sk z j z ■ The next step is to obtain the equations for A^ 
and Bf 1 ' . The mode j is coupled to modes s that satisfy 2wk 



(n) 



\UJj ± UJ S 



Since = 3 Wk this relation is satisfied 
for lo s ~ Wk (as expected) and for lo s — 5oJk- We assume that the spectrum under consideration does not satisfy the 
latter. In this case, the equations read 



dA) 



(n) 



(-1) 



dr 

dB ) n) 
dr 



3 



13* B (n) 



(38) 
(39) 



G 



In order to find the solution to the above equations we write the system in matricial form 



where 



v(t) 



B^(r) 



M = 7 



and the initial condition reads 



-k x 
-k x 
_(_l)*.+i-j s / 3 o 

_(_i)fc»+i^ x / 3 



(40) 










(41) 



v(0) 



fvh=J»A 
) 





1 .5„ 







(42) 



The solution is easily obtained after diagonalizing M. . The eigenvalues are given by 

A = ±^±!V|9A2-12£|, 

where we have used Eq.(^). The solution can be formally written as 

v(t) =Ce T>T C- 1 v(0), 



(43) 



(44) 



where T> is the eigenvalues diagonal matrix and C is the corresponding eigenvectors matrix. This means that A^ ] and 

are linear combinations of exponential functions of the eigenvalues in Eq.(^) times r. The exponential growth 

of and A^ is determined by the eigenvalues with positive real part. Looking at Eq.(Q) we conclude that the 
number of created photons, in both the mode k and the mode j, will increase exponentially in time with a rate given 
by ^k x . 

This is our main result. In the 'resonance parametric case' the resonant mode may be coupled to some other mode. 
In this case the number of created photons in both modes grows exponentially in time with the same rate, which is 
exactly one half of the rate expected for the resonant mode when the coupling is neglected. 

We have derived Eq.(f4f]|) assuming that only two modes are coupled. If the spectrum contains one mode s such that 
lj s = 5wk, besides = 3cjk, we get three coupled modes. The resulting equation will be similar to Eq.(^fj) but with 
a 6 x 6 matrix to diagonalizc. The number of photons in each mode (s, j and k) will grow exponentially in time. Due 
to the non-equidistant character of the spectrum, it is not common to have three modes coupled. For a cubic cavity 
(i.e. L x — L y = L z = L), the first three modes coupled are k = (11, 16, 13), j = (67, 16, 13), and s = (115, 16, 13), the 
frequency of the lowest mode being an order of magnitude larger than the fundamental frequency of the cavity. This 
case is therefore of less interest. 

Let us now discuss briefly what happens with the remaning cases in which the MSA is non trivial and differs from 
the naive perturbation approach. We first study two non- resonant modes s and p (i.e. lu s ^ 0/2 ^ w p ) satisfying the 
condition (E61 



= £1. H. After some algebra on Eqs.(|24|) and (|25|) we get for A^ and Bp 1 ' the following 



? (n) 



dA 



(n) 



dT 



2w s V L 



dB; 



(n) 



dT 



2oj p \ L r , 



(45) 
(46) 



4 Note that we are not necessarily within the 'parametric resonance case', 
frequency of some other unperturbed mode. 



The external frequency Q could be or not twice the 
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and same equations holding for B s and A p . The solutions are straightforwardly obtained, giving for the average 
value of the number operator 

(.V.JH.V^.rff^) 2 ^^ ,47, 

Note that if we set s x — p x (s = p) we recover the 'parametric resonance case'. This example shows the possibility of 
obtaining exponential growth of photons in modes which are not in resonance with the external frequency. 

If the spectrum contains some sequence of equidistant frequencies w Pi separated by SI, the corresponding modes 
will be coupled through the conditions f2 = 0J Pi+1 — LU Pi and Q = uj Pi — w Pj _ r One can show that, as long as the 
modes Pi are not coupled to modes outside the sequence, the number of created photons in each of these modes will 
be an oscillatory function of time. 



B. Examples 



The first and more important example is the cubic cavity. In order to obtain 'parametric resonance' we fix f2 as 
twice the lowest cavity frequency, 



L 



(48) 



For this example we will assume that L — lcm. The fundamental mode k = (1, 1, 1) will be coupled to j = (5, 1, 1) 
because Wfs.^i) = 3^(1^.1). Only these two modes are coupled, since there does not exist in the spectrum any mode s 
satisfaying tv s = 5^(1^.1). The exponential growth for the modes k and j will be one half of that expected by previous 
papers [JL2|. Now we can write explicity Eq.(p4|) for this particular case. The result is 



flf->(r) 



Aj»>(T) 



<5n,k 

\/2^k 



[cos(2.56t) cosh(0.45r) + 0.176 sin(2.57r) sinh(0,45r)] 



[1.76 sin(2.57r) cosh(0.45r)] , 



[cos(2.56t) sinh(0.45r) + 0.176 sin(2.57r) cosh(0.45r) 



[1.76 sin(2.57r) sinh(0..45r) 



-0.586 sin(2.56r) cosh(0.45r)] 

s 

' " ' cos(2.56t) cosh(0.45r) - 0.176 sin(2.56r) sinh(0.45r)], 



/2wT 



-0.586 sin(2.56r) sinh(0.45r)] 
5. 



^=[cos(2.56t) sinh(0.45r) - 0.176 sin(2.56r) cosh(0.45r)]. 



/2lj 



(49) 



(50) 



(51) 



(52) 



An important remark is that this solution satisfies the unitary condition for the Bogoliubov transformation (|ll]) , 

Ei4 n) i 2 -i4 n) i 2 - — 



2w k 



(53) 



We can compute the number of created photons in each mode inserting Eqs.(50) and ([52|) into Eq.(14). The result is 
(Mc) = cos 2 (2.56T / )sinh 2 (0.45r / ) + 1.06 sin 2 (2.56T/) cosh 2 (0.45T/) + 0.088 sin(5.12r/) sinh(0.9r/), (54) 



(A/j) = cos 2 (2.56T/)sinh 2 (0.45r/) + 1.06 sin 2 (2.56T/) cosh 2 (0.45T/) - 0.088 sin(5.12r/) sinh(0.9r/). 



(55) 



When Tf > 1 these expressions are approximated by 
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(Mc) w (TVj) « e°- 9T/ . 



(56) 



In a previous paper |9J the authors considered two dimensional cavities, which means that one of the cavity's 
dimensions is much smaller than the others (say L z << L x , L y ). We can easily recover this limit by omitting the 'z' 
dimension. In what follows we will discuss this case, for increasing external frequencies. 

Let us first assume that L x = L y . If £1 = 2cj( 11 ), then the fundamental mode k = (1, 1) does not couple to any other 
mode and it grows exponentially in time. The next resonant frequencies are £1 = 2uj( 12 ) and £1 — 2cl>( 2 .2)- ln both 
cases the resonant mode is not coupled. If £1 = 2w( 1 3 ), the mode (1, 3) will be coupled to the mode (9, 3) and both will 
grow exponentially in time. However, the mode (3, 1) also satisfies the 'parametric resonance condition' and, being 
uncoupled to other modes, it will grow faster than the previous ones. For the same frequency, £1 = 2w( 1 3), we have 
found by inspection three equidistant modes, ^(13,39) = 13ny/lO/L, 1^(27,39) = 15tt>/T0/L and 0^37,39) = 17iry/l0/ L. 
The number of created photons in each mode will oscillate in time. 

Now we choose L x = 3L y . This choice makes the fundamental mode to be coupled in 'parametric resonance'. If we 
set f2 = 2w(n), the mode (9, 1) satisfies that W(g,i) = 3w( la ), so both modes will grow exponentially. 



C. Off resonance 



In this subsection we study what happens when the external frequency £1 is slightly off-resonance, i.e., O = fi + h, 
where £1 is a resonant frequency and h <C £1. We will assume that h = ea, where a = 0(£1). We will show how to 
apply MSA to this case. Off-resonance motions have already been considered in the literature [jl8] using a different 
approach, and it was shown that there are threshold conditions on h for exponential photon creation. 

For small amplitudes of oscillations, the equations for the modes are still Eq.(16) with £1 replaced by the external 
frequency f2. Since h <C £^, all factors of the form e ±lht may be regarded as slow oscillations, so that the MSA 
conditions to get rid of secula rities are exactly the same as in the resonant case, Eqs.( prjj27| , p^ ). However, the 
equations for the modes 



do get modified. Their off-resonant version reads 



dA 



(n) 



(It 



2 Wk L2 

E 



e*°" B^S(2^ £1) + + ^±h.)S{- U]l - Wj + £l)^ 5kje — 



tt + h 



2w k 
j 

— — )6(wk - wj - n) + e MT (wj — )S(uj k - «j + fl) 



n + h 



(57) 



and 



dB 



(n) 



dr 



E 



4 n) ^(2 Wk -£l)+^(- 



Q + h 



2 )*("«*- "J + «)- 2£|/k 



(n) 



2w k 



(n) 



(58) 



Let us solve these equations in the quasi parametric resonant case, that is for Q — h = $1 = 2w k . In the case when 
there are no coupled modes, we get two coupled first order differential equations for and B^\ After the change 
of variables 



.(n) _ iar/2 (n) 
~~ c "k > 

r>( n ) _ iar/2i(n) 
D k — e u k ) 



(59) 



the equations take the form 



(n) 



dr 



db 



(n) 



(n) 



* 2 k 



7fe x a k . 



(60) 
(61) 



There will be growing exponential solutions if an eigenvalue of the corresponding matrix 
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has a positive real part. The eigenvalues are given by 



and lead to the following threshold for resonant behavior 



\h\ „ e 



(fc 



< 2e 7 k x <=> y < (64 ) 



(fc)' 



For this uncoupled, quasi-resonant situation, the resonant mode satisfies a Mathieu equation (see Eq.(|16[)). The 
threshold we obtained in Eq. ( |6~4| ) coincides with the one obtained for the Mathieu equation using a different method 
& 

As a second example, let us consider the case of two coupled modes, say modes k and j, for which lo- } = 3o>k- After 
the change of variables ([59|), and defining 



(n) _ 3iar/2 (n) 
j ~~ e "j ' 

( n ) _ -3iar/2i (n) 



B w = e -s.or/a & W i (65) 
the off-resonance form of Eqs . (|35| J3^ , ^8| , ^9|) is 

^ - «L n) + 7 (-fc, &L n) + (-iy-+*-jx «i n) ) , (66) 



(n) 



+ 7(-fc,a^ n) + (-l)^+Sx&j n) ), (67) 



j_ = _3i| a f - (-l)^f aW, (68) 

^ = 3z|6f)-(-l)^^^, (69) 

where we have neglected terms proportional to h/uu^ = 0(e) in the r.h.s., since they would introduce only small 
corrections to the eigenvalues and thresholds. Just as in subsection A, we can rewrite these equations in matricial 
form, and find the solution by diagonalizing such a matrix. The corresponding eigenvalues are 



where U = + ^[2^ — 3fc^] is a positive number (see Eq. (|37|)), and 

V = 4a 4 + Aa^(kl + \f x ) + 7 4 fc^2 _ ^ 

For a = 0, V is negative, so there are two eigenvalues with positive real part, which corresponds to exponential 
growth. As \a\ grows, V will eventually become positive, all eigenvalues are pure imaginary for \/V < U, and no 
exponential growth is obtained. The condition for exponential growth is V < 0, which again sets a threshold for the 
off- resonance frequency difference h, 



\h\ <e y\l-di + x/4 + d 2 , (72) 
with 

di = (k 2 x + > 0, (73) 

da = A£(~£ - Ag) > 0. (74) 

For the example of the cubic cavity discussed in the previous subsection, the threshold is \h\ < 0.68e7, that is 
\h\/n < 0.06e. 
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V. DYNAMICAL CASIMIR EFFECT AT FINITE TEMPERATURE 



Up to this point we have assumed that the field was in the in- vacuum state T = (see, for example, Eg. (|l4|) ) . 
ft is well known that the temperature contribution can dominate the pure vacuum effect when computing the static 
Casimir force |pc)| . Then, it is expected that temperature plays an important role in the dynamical Casimir effect as 
well. In what follows we shall derive an expression for the number of created photons inside the cavity equivalent to 
Eq. (|l4|) , but now assuming an initial state in equilibrium at finite temperature. After that, we will apply the obtained 
result to the resonant vibrating cavity. 

We adopt the scheme of scalar field quantization developed in section II. For t < the wall is at rest, so we assume 
the system to be at thermal equilibrium at finite temperature T = 1/(3. The state of the system is described by a 
statistical operator p which does not evolve in time. We expand this operator in Fock states of the field at t < 

P=\ E E ••■e- /3i: - ( " k - + ^ )Bl -hk 1 n k2 ...)(n kl n k2 ...|, (75) 

where = lu^ and the normalization factor is given by 

Z = E ...e-^i^+i^. (76) 

nk x >0 nk 2 >0 

At t = the wall begins to oscillate and the system will no longer remain at thermal equilibrium. Following the 
steps given in section II we assume that the wall stops at t = ifi na i- In that case, the number of photons in the k-mode 
is the average value of a^'a^"* with respect to the initial state p 

(Mc)p = Tr(pat° ut ar). (77) 
Using Eq. ( [ll| ) we can write this expression in terms of a ln and 111 

Tr(p o£ out ak *) = ^2 /?„ k aj k Tr(p a^af) + /3„ k /3f k Tr(p a^aj m ) 

nj 

+a£k «jk Tr(p&i ia a?) + a* k /3f k Trfp^a]" 1 ). (78) 
With the help of Eqs.(^5|) and ( |76|) it is a straightforward calculation to find that 

Tr(p a*af) = = Tr(p a^aj in ) (79) 
Tr(p al ia af) = _ - 5 ni = Tr(p a„ n a] ,n ) - 1. (80) 



Inserting this into Eqs.(|7q) and (77) we arrive at 

<A/k) P - £ l/^l 2 + Ed^kl 2 + |a nk | 2 ) ^r— r (81) 

ii n 

Finally, using Eqs.([L3|) and ( |l4| ) we get 

(■Mc)p = (AA k )T=o + Ed^l 2 + l^l 2 )^^- (82) 

Let us now apply the results obtained in section IV for the vibrating cavity to the case in which the field is initially 
at thermal equilibrium. In the parametric resonance case without coupling the Bogoliubov coefficients are diagonal 
(see Eqs.(|2|) and @)), so Eq.(|^) can be reduced to 

(-Mc)fl - (A/k)T=o (l + 2 * - ) + * - , (83) 



with (A/"k)T=o given by Eq.(|34|). We see that the effect of the temperature is to enhance the amount of created photons 
in the pure vacuum case by a thermal distribution factor. Note that the second term in expression (83) corresponds 
to the average number of photons, in the mode k, present in the cavity when t < 0. The same result was obtained in 
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Ref. |L2| for the fundamental mode of a cubic cavity using a different approach. However, we have seen that in this 
case the Bogoliubov coefficients are not diagonal, due to the coupling between the fundamental mode k = (1, 1, 1) 
and the mode j = (5, 1, 1). With the help of Eqs.(fi"9|), (50), (|5l"|), and ( |5^ ) it is easy to obtain the number of photons 
present in the fundamental mode after the wall stops 

(7Vk)p = (Mc)r=o(l + 2n k ) - 2sin 2 (2.56r) cosh 2 (0.45T)(n k - tij) 

-sin 2 (2,56r)(n k + nj) + n k , (84) 

where n k = l/(eP Ek — 1) is the Bose mean ocupation number for t < 0. For the mode j we find 

(Mj) p = (JVj) T =o(l + 2nj) - 2sin 2 (2.56r)cosh 2 (0.45r)(nj - n k ) 

-sin 2 (2.56r)(nj + n k ) +rzj. (85) 

Again, we obtain that the effect of the temperature is to increase the pure vacuum case with thermal factors now 
depending on t. For Tf > 1 we have (A/k)r=o ~ (A^)t=o ( see Eq. [36j). Therefore the total number of photons created 
inside the cavity becomes 

(Motal)p = (1 + ™k + nj)(Motal>T=0- (86) 

For a cubic cavity of size L — 1cm at room temperature T w 290K, we obtain that the total number of created 
photons is approximately 300 times that of pure vacuum (T = 0). 



VI. CONCLUSIONS 



We have calculated the photon production inside a three dimensional oscillating cavity, using MSA to deal with 
the resonant effects. We have taken into account that, even though the spectrum of the cavity is non equidistant, the 
different modes may be coupled, and this coupling affects the rate of photon creation. 

We have found resonant effects when the external frequency is equal to the sum of the frequencies of two unperturbed 
modes = lu s + uj p . When lu s =/= lo p , the number of photons in both modes grows exponentially. When Q = 2o> k , 
the usual "parametric resonance case" , the number of photons in the mode k also grows exponentially, along with 
the number of photons of other modes coupled to k. When the mode k is coupled to one mode, the rate of photon 
creation decreases by a factor 2 with respect to the uncoupled case. 

We have also analyzed slightly off-resonance situations. Using an extension of the MSA we showed that the number 
of photons in the relevant modes also increases exponentially if certain threshold conditions are satisfied. These 
conditions imply that the external frequency should be almost equal to the resonant frequency. 

As an important example, we have described in detail the case of a cubic cavity. The fundamental mode (1, 1, 1) is 
coupled to the mode (5, 1, 1) when the external frequency is Q = 2Lj^ i iy The number of photons created in both 
modes grows as e°' 9T . Neglecting the mode coupling, one would erroneously conclude that the number of photons in 
the mode (1, 1, 1) grows as e 18T , and that there is no exponential growth for the mode (5, 1, 1). The mode coupling 
in three dimensional cavities has not been taken into account in the previous literature. 

We have also computed the enhancement of the dynamical Casimir effect for an initial thermal state. The main 
result is contained in Eq. (|82|). Only when the Bogoliubov coefficients are approximately diagonal one recovers the 
usual result, i.e. the number of photons created in a given mode at temperature T is the T = result times the 
thermal distribution factor. 

For simplicity, we have considered a quantum scalar field. The generalization to the case of an electromagnetic 
field is not completely straightforward. We state here the main results, the details will be described in a future work. 
Assuming that the potential vector satisfies the gauge condition V.A = 0, the Maxwell equations read DA = 0. The 
boundary conditions are the usual (perfect conductors) on the static walls. On the moving mirror, these boundary 
conditions must be imposed in a Lorentz frame in which the mirror is instantaneously at rest (see |lC| ) . This implies 
that, at x = L x (t), the potential vector must satisfy A y (t, L x (t), y, z) — A z (t, L x (t),y, z) = 0. 

As for the case of the scalar field, one can expand the potential vector in an instantaneous basis, and use MSA to 
deal with the secular terms. The resonant conditions coincide with Eqs. ^ - 12^ in the present paper. The differential 
equations for the (now polarization dependent) Bogoliubov coefficients are different for the transverse electric TE and 
transverse magnetic TM cases. However, it can be shown that, in the parametric resonant case, the rate of photon 
production in both polarizations TE and TM is the same and coincides with the rate for the scalar field computed in 
this paper. 
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Finally, we would like to comment briefly about the possibility of observing experimentally the dynamical Casimir 
effect. Considering cavities of dimensions of the order of 1cm, the external frequency should be at least 1GHz in 
order to have resonant photon creation. This is not trivial, the upper limit being hopefully around 100MHz with the 
present techniques |^T| . Other serious technical problem is that, as already mentioned, in order to have resonance the 
external frequency should be tuned with high accuracy to the resonant frequency. Although extremely difficult, an 
experimental verification of the photon production seems not completely unrealistic. 
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